(N 

o 

(N 



C/3 



Families of superhard crystalline carbon allotropes induced via cold-compressed 
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We report a general scheme to systematically construct two classes of structural families of super- 
hard sp^ carbon allotropes of cold compressed graphite through the topological analysis of odd 5+7 
or even 4+8 membered carbon rings stemmed from the stacking of zigzag and armchair chains. Our 
results show that the previously proposed M, bct-C4, W and Z allotropes belong to our currently 
proposed families and that depending on the topological arrangement of the native carbon rings nu- 
merous other members are found that can help us understand the structural phase transformation 
of cold-compressed graphite and carbon nanotubes (CNTs). In particular, we predict the existence 
of two simple allotropes, R- and P-carbon, which match well the experimental X-ray diffraction pat- 
terns of cold-compressed graphite and CNTs, respectively, display a transparent wide-gap insulator 
ground state and possess a large Vickers hardness comparable to diamond. 

PACS numbers: 61.50.Ks, 61.48.De, 64.60.My 



a 



o 



> 

oo 

rn 

O 



X 



Carbon exhibits numerous allotropes (fullerenes, car- 
bon nanotubes(CNTs), graphene, graphite, diamond and 
amorphous carbon) thanks to its ability to form sp-, 
sp^- and sp^-hybridized bonds [lH3|- It is well-known 
that compression of graphite at high pressure (>15 GPa) 
and high temperature (> 1300 K)|^] leads to the forma- 
tion of cubic or hexagonal diamond. In contrast, cold 
compression of graphite [4l4l5l| , single- walled and multi- 
walled CNTs [16i-t22i] results in superhard allotropes of 
carbon, which were found to be intrinsically different 
from hexagonal (or cubic) diamond |5l. Il6|. Upon pres- 
sure release, the obtained cold compressed graphite phase 
can be quenched only at low temperature (< 100 K) [1] 
whereas CNT phases compressed above 75 GPa can be 
quenched at room temperature [l6|. In addition, these 
high pressure phases exhibit superior mechanical perfor- 
mance with the ability to indent single-crystal diamond 
[a, lla] , indicating at least comparable hardness to dia- 
mond. 

However, the experimental crystal structures of these 
cold compressed phases of graphite and CNTs remain 
heavily debated. In an effort to shed light on the puz- 
zling structural problem of cold compressed graphite sev- 
eral superhard sp'^-hybridized candidates (monoclinic M 
carbon [23|, |2J|, body-centered tetragonal bct-C4 carbon 
illii], orthorhombic W [27J and Z carbon Hll^l) have 
been proposed. To date, among all these proposed al- 
lotropes the most stable one is Z carbon [29], computa- 
tionally predicted through the minima hopping method 
(MHM), which is exactly the same as the 0CI6-II phase 
proposed recently by metadynamics simulations of struc- 
tural transformations [28| and the C-centered orthorhom- 
bic Cg phase proposed by particle swarm optimization 



(CALYPSO) [20] • Z carbon was thought to be the best 
candidate [23] so far that can explain the experimental 
XRD peaks and Raman active mode for cold compressed 
graphite [a, l29| . Although Zhao et al. [30| argued there 
was a matching problem with the XRD of cold com- 
pressed graphite, they claimed [30(| that the Z (namely, 
Cg) phase can be interpreted as the structural solution of 
the quenchable superhard carbon phase recovered from 
cold compressed CNT bundles [16|. However, we found 
that its XRD patterns still differ significantly compared 
to the experimental data (as discussed below). 

If we pay more attention to the previously proposed al- 
lotropes, both W and M carbon phases can be described 
as corrugated graphite sheets interconnected by an al- 
ternating sequence of odd 5+7 membered rings of car- 
bon 28|. Similarly, both bct-C4 and Z carbon can be 
characterized by alternating even 4+8 membered rings. 
Considering all possible even and odd rings of carbon 
(4+8 membered rings in bct-C4 and Z carbon, 5+7 mem- 
bered rings in W and M carbon, 6-membered rings in 
diamond), many additional possible combinations could 
be expected. Therefore, it would be highly desirable to 
seek a general scheme to understand the principles of the 
structural formation of possible alternative superhard al- 
lotropes related to cold compressed graphite and CNTs. 

In this pa per , through first-principles calculations (for 
details, see [31]), we report on two S and B families of 
sp^-hybridized superhard carbon allotropes by discussing 
the topological stacking of zigzag and armchair carbon 
chains consisting of odd 5+7 and even 4+8 membered 
ring patterns (see Fig. SI [3l|). Our analysis demon- 
strates that, after introducing the hexagonal rings which 
separate the periodic 5+7 or 4+8 membered patterns, a 



series of new structures can be readily created. We fur- 
ther elucidate the energetic, mechanical and electronic 
properties of the obtained novel phases, confirming that 
all these phases are wide-gap transparent and superhard 
insulators. In particular, our currently proposed P car- 
bon as well as several other new phases are energetically 
more favorable than the previously known most stable 
Z carbon 28|, |30| over a large range of pressures. More- 
over, we confirm that R carbon of the S-family and P car- 
bon of the B-family are the most likely candidates so far 
to match the experimental data of the cold-compressed 
graphite [5| and CNTs [16|, respectively. 

Each graphite sheet (i.e., graphene) can be described 
as an ordered array consisting of infinite long zigzag 
or armchair chains of carbon atoms. The direct con- 
nection between zigzag and armchair chains results in 
a tilt grain boundary [32| which is interconnected by 
5-membered and 7-membered rings of carbon, different 
form the ideal 6-membered ring in graphene. Similarly, 
the two-dimensional (2D) projections along the c-axis for 
M carbon and W carbon of cold-compressed graphite [29| 
are also composed of 5-1-7 ring patterns, realized by the 
repeated stacking of zigzag and armchair chains of carbon 
in Fig. S2 3l|. The most important feature of both M 




and W carbon is that each armchair chain or each zigzag 
chain is not infinitely long. It is recognized that carbon 
bonds are combined in group of four armchair or zigzag 
bonded chains. As illustrated in Fig. S2 [3l|, these four- 
bonds chains are interconnected by a step. Interestingly, 
if the directly connected zigzag and armchair chains in M 
carbon would be infinitely long (without any step), the 
simple S carbon phase could be easily realized. As shown 
in Fig. la, its 2D-projection can be described as the re- 
peating parallel array of directly connected zigzag and 
armchair chains, consisting of diagonally opposed 5-1-7 
ring pattern. If this 5-1-7 ring pattern is further twofold 
rotated, another new structure called R carbon can be 
formed (Fig. lb). Therefore, the unit cell of R carbon is 
twice the size of S carbon. 

The structures of S carbon and R carbon can be fur- 
ther modified. As shown in Fig. Ic and Id, by inserting 
two zigzag or armchair chains to separate the 5+7 ring 
pattern in S carbon, two new structures (S-S1Z2 and S- 
S1A2) can be realized. Here, the notation of S come from 
the 5-1-7 ring patterns in S carbon (Si means only one S 
unit), whereas Z and A denote the infinitely long zigzag 
and armchair chains (Z2 and A2 refers to two zigzag and 
armchair chains), respectively. Furthermore, it is possi- 
ble to continuously extend the number of Z and A chains 
to construct a series of new structures named S-SmZ2„ 
and S-SmA2„. Here, the parameters m and n are integers 
(0, 1, 2 •••). Following this nomenclature, R carbon can 
be described as S-S2 in which S' denotes the existence of 
the twofold rotated symmetry of 5+7 ring pattern at vari- 
ance from the normal 5+7 ring in S carbon. Following 
similar considerations, S' can be further extended into a 
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FIG. 1: (color online) (a), (b), (c) and (d) illustrate the crys- 
tallographic structures of the proposed S carbon, R carbon, 
S-S1Z2, and S-S1A2 phases. Left column: three-dimensional 
crystalline structures and right column: the two-dimensional 
(2D) projections of the 2x2x1 supercell along their b- or c- 
axis. For S carbon, the 3x1x2 supercell is used. The grey 
background highlights the armchair chains and white back- 
ground denotes the zigzag chains. 



FIG. 2: (color online) (a), (b) and (c) illustrate the crystal- 
lographic structures of the bct-C4 carbon, Z carbon, P car- 
bon (this work). Left column: three-dimensional unit cells 
and right column: the 2D projections of the 2x2x1 super- 
cell of bct-C4 and of the 1x2x1 supercells of Z carbon along 
their c-axis, and the 2D projections of the 1x1x2 supercells 
of P carbon along its b-axis. (d), (e) and (f) show the 2D- 
projections in the 1x2x1 supercells along the c-axis for three 
typical lattice structures of the B-B1A4 series. 



TABLE I: DFT optimized lattice constants (a, 6, and c in A) 
Wyckoff position(W.p.), bulk and shear moduli (B and G in 
GPa), estimated Vickers hardness (H„ in GPa) and theoreti- 



cal density (/ 


J, g/cm ) for S, R and P allotropes. 


Types 


W.p. X y z 



S carbon a 

(P2/m) b ■ 
c - 

P 



4.7302 2m 

2.4950 2m 

4.0837 2n 

106.1° 2n 



0.1175 
0.5344 
0.1131 
0.4209 



0.0 
0.0 
0.5 
0.5 



0.6746 G = 
0.3333 B = 
0.8997 H^ 
0.1319 p = 



457.7 

412.6 

= 78.3 

3.44 



R carbon a - 


= 7.7886 4g 


0.6731 0.9630 


0.0 


G = 462.4 


(Pbam) b = 


= 4.7752 4g 


0.8435 0.8087 


0.0 


B = 434.2 


c - 


= 2.4958 4h 


0.9546 0.8613 


0.5 


H^ = 75.0 




4/1 


0.5704 0.8926 


0.5 


p = 3.45 



P carbon a 

{Pmmn) b 

c 



8.6650 4/ 

2.4875 4/ 

4.2160 4/ 

4/ 



0.5357 
0.2077 
0.0414 
0.7151 



0.25 
0.25 
0.25 
0.25 



0.4322 G = 

0.4348 B = 

0.5625 Hy 

0.4343 p = 



485.0 

449.1 

= 78.5 

3.51 



series of new structures called S-Sj„Z2n and S-S^A2n. It 
should be noted, that when n is increased to n+l the 
number of carbon atoms in the unit cell is increased by 
eight atoms, due to the fact that eight carbon atoms can 
form a full cycle with a sixfold ring in the projections 
along the other b or c-axis. As such, the atom numbers 
of the conventional unit cells of these families are equal 
to 8{m + n). 

Figure 2a shows the 2D projection of bct-C4, indicat- 
ing the simplest 4+8 rings pattern along the armchair 
orientation (direct face-to- face opposed connection of two 
armchair chains). By inserting the armchair chains into 
periodically separated 4+8 ring pattern, a new structural 
family of B-BmA2„ can be constructed. Here, B is the 
unit of the 4+8 rings composed of two opposed armchair 
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FIG. 3; (color online) Calculated enthalpy difference per atom 
with respect to graphite for a series of carbon allotropes as a 
function of pressure. The stable pressures of these proposed 
allotropes compared to graphite are labeled. 



chains originated from bct-C4 and A refers to the arm- 
chair chain. For the sake of convenient comparisons, we 
further split A2„ into h-i^nijin" to distinguish the left- 
and right-armchair chain numbers of the B unit. Clearly, 
n'+n" is equal to 2n, indicating that the number of atoms 
is 8(m+n) in the conventional unit cells. With m—1 and 
n =0, B-Bi becomes bct-C4. Furthermore, the addition 
of two armchair chains with n=l separates the 4+8 ring 
pattern thus leading to two new structures: the first one 
(Fig. 2b) corresponds to the Z carbon (B-BiAliai) as 



proposed in Refs. |28l - l30l |. whereas we name the second 
one (Fig. 2c) P carbon (B-BiALoiJ2)- By adding four 
armchair chains, n=2, three alternative structures are 
formed, as depicted in Fig. 2d (B-BiAi2fi2), Fig- 2e 
(B-BiAlij^3) and Fig. 2/ (B-BiAiofl4). Therefore, it is 
clear that all structures in the B-B„iA2„ family are com- 
posed of 4+8+6 ring patterns. The introduction of zigzag 
chains in the 4+8 ring patterns bct-C4 will lead to the oc- 
currence of 4+8+5+7 ring complex structures. However, 
we do not discuss these cases since all phases associated 
with 4+8+5+7 ring complex (not shown here) are ener- 
getically much less stable than the S- and B-families. 

Figure 3 compares the relative thermodynamic stabil- 
ities of all proposed allotropes. We find that, for both S 
and B families with increasing zigzag or armchair units 
(namely, increasing the 6-membered rings in the unit cell) 
their structures become more stable in energy and the 
transition pressure associated with the structural transi- 
tion decreases. This can be attributed to the introduction 
of more ideal 6-membered rings which actually reduces 
the total strain in system. If the 5+7 patterns (or 4+8 
patterns) are doubly coupled in a complementary mode, 
the energy of system can be even lower. The complemen- 
tary mode reduces the strain of system compared to the 
cases in which a single 5+7 (or 4+8) pattern appears in 
a non-coupled complementary one. Although Z and P 
phases share the same composition of B-B1A2, P carbon 
is more stable in energy than Z carbon (Fig. [3]) because 
the 4+8 patterns in P carbon are always coupled doubly 
in a complementary mode (Fig. [2t) , whereas in Z carbon 
they are separated by the 6-membered rings (Fig. Wp)- 
In addition, in these families as n increases the energies 
of the structures approach more closely that of diamond 
due to the increase in the number of ideal 6-membered 
rings. In this situation, the S (5+7 ring) and B (4+8 ring) 
patterns can be thus viewed as the defect or tilt grain 
boundary in diamond. Following this viewpoint, the S- 
SmZ2n and S-SmA2„ families can be interpreted as the 
combination of diamond (6-membered ring) and S carbon 
whereas the B-BmA2n families show the combined char- 
acter of diamond and bct-C4. Furthermore, from Fig. [3] 
our selected phases P carbon, B-BiALoi?,4, B-BiA^ifla, S- 
S1Z2, and S-S1Z4, are energetically more favorable than 
all previously theoretically proposed structures includ- 
ing Z carbon. In particular, even when their vibrational 
entropies and zero-point energies derived from phonon 



densities of states are taken into account, the relative 
stabilities for S-, R- and P-carbon allotropes remain un- 
changed at least in the temperature range from to 800 
K at both and 15 GPa. 

Structural optimization revealed that R and P carbon 
crystallize in the orthorhombic structure whereas S car- 
bon has a monoclinic structure (Figs. 1 and 2). The 
structural details of S, R, P and the other nine new al- 
lotropes proposed here are listed in Tab. 1 and Tab. SI 
3lj , respectively. The phonon dispersions in the whole 
Brillouin zone for S, R and P carbon phases have been 
derived, confirming their crystalline stabilities even at 
GPa (Fig. S3 31]). In addition, our calculations uncov- 
ered that all these allotropes exhibit large bulk and shear 
moduli (Table 1) which are comparable to those of dia- 
mond. Utilizing our recently proposed formula 33|, |34 j , 
the Vickers hardness of S, R and P carbon are estimated 
to be 78.3, 75 and 78.5 GPa in Tab. 1, respectively (de- 
tails refer to Tab. S3 ^31]. This clearly demonstrates 
that these phases are superhard in agreement with the 
quasi- sp'^-hybrid covalent bonding framework. The HSE 



electronic band structures illustrate that S-, R-, and P- 
carbon are all wide-gap transparent insulators with band 
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FIG. 4: (color online) (a) Simulated XRD patterns (with 
wavelength of 0.3329 A) for selected phases at 18.4 GPa 
compared with experimental results from Ref. [a| for cold- 
compressed graphite, (b) Simulated XRD patterns of the five 
allotropes in B family at ambient pressure, in comparison with 
the date measured from the quenchable phase recovered from 
cold-compressed CNTs tlQ]. The red arrows highlight five 
main experimental peaks. 



gaps ranging between 3.9 and 5.5 eV (Fig. S3 [3l|). In 
particular, both S and R carbon have the widest direct 
band gap at F (about 5.5 eV) among all known allotropes 
discussed here. 

Although Z carbon was identified as a good candidate 
of cold compressed graphite in the latest work [29[ , this 
was challenged by Zhao et al. [30| due to the prob- 
lems for matching experimental XRD patterns. They 
demonstrated that its simulated patterns match neither 
the experimental peak at 16.8° and nor the gradually de- 
creased peak density around 9.0° [30| for cold-compressed 
graphite (for details, see the online supporting materials 
of Ref. [29[). Importantly, we found that the simulated 
XRD patterns of R carbon match much better the exper- 
imental data for five main peaks as marked by arrows in 
Fig. mj upper panel). The obtained pressure of 11.5 GPa 
for the stability of R carbon is consistent with the exper- 
imental value. Therefore, R carbon can be considered as 
a highly likely candidate for cold-compressed graphite. 

Furthermore, we found that P carbon is a better candi- 
date to interpret the experimental XRD pattern of cold- 
compressed CNTs than Z carbon as claimed in Ref. [30[ . 
From Fig. 2] (lower panel) , for Z carbon and three other 
B-family phases (B-BiAL2i?2, B-BiA^ii^a, B-BiALofl4) 
the simulated peaks at d — 2.155 A are stronger (or at 
least comparable) than the ones at 2.053 A. This no- 
tion is not consistent with the reported XRD patterns 
|16| which clearly revealed the d=2.053 A peak has the 
highest intensity. This fact suggests that Z carbon may 
not be a good candidate of cold-compressed GNTs. Nev- 
ertheless, it is exciting to note that the simulated XRD 
pattern of P carbon shows a significantly improved agree- 
ment with the experimental data. Its five main peaks 
at d = 2.167 A, 2.107 A, 1.510 A, 1.244 A and 1.179 
A give satisfactory accordance in both locations and in- 
tensities with experimental d spacings of 2.155 A , 2.053 
A, 1.495 A, 1.248 A and 1.161 A respectively (see Fig. 
m^lower panel)). Besides, its theoretical bulk density (p= 
3.51 g/cm'^) and bulk modulus {B = 449.1 GPa) are in 
nice agreement with the experimental data {p = 3.6±0.2 
g/cm3 and B = 447 GPa), see Tab. 1 and S2 [sij. 

In summary, R-carbon in the S-family and P-carbon in 
the B-family have been found to best match with the ex- 
perimental data of cold-compressed graphite and CNTs, 
respectively. Moreover, for all other proposed allotropes 
related to the S- and B-families, the simulated XRD pat- 
terns also capture their main experimental features. In 
general, with increasing numbers of zigzag or armchair 
units, the main peaks remain essentially unchanged (see 
Fig. S4 3l|), though their corresponding 29 values are 
slightly shifted. Considering that in real samples a pres- 
sure gradient could exist [a,[la| and the tube diameters of 
CNTs show a wide range from 1.8 to 5.1 nm [16] which 
may create the conditions of the formation of different 



phases 25|, |35J, it would be reasonable to expect that 
the cold-compressed phase of graphite and CNTs can be 
interpreted as a mixture of several of the proposed al- 
lotropcs. 
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